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ABSTRACT 

We present a deep 1.4 GHz survey made with the Austraha Telescope Compact Ar- 
ray (ATCA), having a background RMS of 9 /xJy near the image phase centre, up to 
25 /j.Jy at the edge of a 50' field of view. Over 770 radio sources brighter than 45 /iJy 
have been catalogued in the field. The differential source counts in the deep field pro- 
vide tentative support for the growing evidence that the microjansky radio population 
exhibits significantly higher clustering than found at higher flux density cutoffs. The 
optical identification rate on CCD images is approximately 50% to R = 22.5, and the 
optical counterparts of the faintest radio sources appear to be mainly single galaxies 
close to this optical magnitude limit. 

Subject headings: surveys — radio continuum: galaxies — galaxies: general 
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1. Introduction 

Our understanding of the faintest radio sources has 
advanced over the last decade through increasingly 
sensitive radio surveys and through follow-up work at 
optical and infrared wavelengths. The appearance of 
a 'new' population of radio galaxies at millijansky and 
sub-millijansky levels, first revealed by the changing 
slope of the 1.4 GHz radio source count distribution, 
has been confirmed by spectroscopic and multicolour 
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il. 1998). The population comprises star forming 
galaxies whose numbers increasingly dominate classi- 
cal (AGN-powered) radio sources as the flux density 
limit falls below 1 mJy. The evolution of these galax- 
ies, and the role of interactions and mergers in the 
population, are only partially understood. 

The number of catalogued microjansky radio sources 
at 5 and 8 GHz is a few hundred, and at 1.4 GHz there 
are few catalogued sources fainter than 100 ^Jy. The 
Phoenix Deep Survey (PDS) aims to increase the num- 
ber of known 1.4 GHz sources fainter than 100/LtJy. 
The PDS covers a 2° diameter field selected to lie in a 
region of low optical obscuration and devoid of bright 



radio sources (Hopkins et al. 1998). The PDS pro- 
vides a large, homogeneous sample of 1.4 GHz sources 
which, through multiwavelength observations, is be- 
ing used to investigate starburst and post-starburst 
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. We present here the source counts and op- 



tical identifications in a Phoenix mosaic image which 
contains many sources fainter than 100 /zJy. 

2. Observations, data reduction and imaging 



Earlier radio observations of the PDS at 1.4 GHz 
provided images with a 5cr sensitivity of 300 /iJy over 
the 2° Phoenix Deep Field (PDF), and 100 /xJy in a 
sub-region of 36' diameter, refer red to as the Phoeni x 
Deep Field Sub-region (PDFS) ([Hopkins et al. 1998|) . 
Additional observations were made in 1997 November 
and December, using the ATCA-6C configuration in 
mosaic mode with 32 channel, 128 MHz bands centred 
at 1.380 GHz and 1.472 GHz. The mosaic has seven 
pointings, six arranged hexagonally around a point- 
ing centred on the PDFS [a(J2000) = 01^ 11™ 13?0, 
(5(J2000) = -45° 45' O'.'O]. A integration of 122 hours 



was spent on the mosaic. The primary calibrator was 
B1934-638. The phase calibrator, B0153-410, was ob- 
served for ten minutes every hour. An additional 42 
hours of data obtained in 1994 September on the cen- 
tral PDFS pointing were included in the imaging. 

The MIRIAD (Multichannel Image Reconstruction, 
Image Analysis and Display) data analysis package 
was used to edit and calibrate the data, as with pre- 
vious PDS observations (Hopkins et al. 1998). An ar- 
tifact at the phase centre of earlier PDS images due to 
self-interference from the sampler clocks is no longer 
present. 

An imaging cell of 2 x 2 arcsec^ per pixel was used, 
providing at least 3 pixels across the shortest FWHM 
of the synthesised beam. Naturally and uniformly 
weighted images were produced. The RMS noise in 
the naturally weighted image is lower, but the synthe- 
sised beam larger (a x (5 « 6" x 12" compared with 
5" X 7" for uniform weighting). Sources have been de- 
tected in the naturally weighted image, and individual 
sources (when detected) have been investigated using 
the uniformly weighted image. 

Mosaic components at each pointing were made 
2100 pixels (70') square, to allow correct CLEANing 
of the sidelobes of out-of-field sources. Although the 
PDF and PDFS fields were chosen to avoid bright 
sources, the two westernmost pointings were affected 
by the sidelobes of a strong source outside the field. 
These were removed by imaging, modelling and re- 
moving this source (using the MIRIAD task UVMOD) 
when constructing the images for those pointings. 
The images for each pointing were CLEANed individ- 
ually, in each case reaching the thermal noise, before 
being combined to form the final mosaic. The central 
50' diameter of the mosaic is shown in Figure |^. 

Because the final image is a mosaic in which the 
individual pointings have different integration times, 
spatial overlap, and primary beam attenuation, the 
distribution of noise in the image is not uniform. The 
RMS noise is a minimum at the centre (9 /iJy) ris- 
ing uniformly to 25 /iJy at a radius of 25'. The RMS 
noise at the centre is a factor of 2.5 lower than the 



previous PDFS image (Hopkins et al. 1998), consis- 
tent with the 4-fold increase in integration time and 
the absence of the phase-centre artifact present in the 
earlier observations. The dynamic range of the image 
can be described in two ways: the brightest source 
(in the NE, 22.7 mJy) lies more than a factor of 100 
above the RMS noise of its surrounds, while the ra- 
tio of the flux densities of the brightest and faintest 
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reliably detected sources exceeds 500:1. 

3. Radio and optical sources 
3.1. 1.4 GHz sources 

The mosaic achieves wide-field coverage and good 
sensitivity at the price of having an unavoidably non- 
uniform noise distribution. The detection and statis- 
tical characterisation of sources is thus a more com- 
plex procedure than it is in single-pointing interferom- 
eter images. An important auxiliary in this process is 
the noise image, in which each pixel is assigned a value 
equal to the theoretical la noise level accounting for 
the observing time, mosaic overlap and primary beam 



atten uation. The miriad task SFIND ([Hopkins et al 



1998 ) was used to compile a list of sources which have 
a peak flux density greater than 4 times the value in 
the noise image at the same position, lying in the cen- 
tral 50' diameter disk of the mosaic. This disk is the 
area bounded by the 25/iJy (la) contour. Over 770 
sources within this region having a peak flux den- 
sity brighter than the local value of 4(t survive vi- 
sual inspection in the interactive confirmation phase 
of SFIND. There is some deviation from Gaussian 
statistics in the distribution of negative pixel values 
0.1% of pixels lie in a non-Gaussian tail associ- 
ated with visible sidelobes and the edges of the im- 
age) which evidently exists in the positive pixel values 
as well, albeit masked by true sources. However, the 
procedure we have adopted provides a very robust 
estimate of the noise level above which we claim de- 
tections, and excludes essentially all spurious positive 
candidates. 

The construction of the source count distribution 



takes account of two significant corrections (c.f. Bop- 
l^ins et al. 1998| ). First, the use of a peak detection 
algorithm to locate sources will underestimate the to- 
tal number of sources to a given total flux density. 
A correction based on the statistical distribution of 
source sizes has been applied. Secondly, it is neces- 
sary to correct the raw counts at any noise level for 
the fraction of the total area over which such sources 
could be detected above 4ct. Figure § illustrates the 
magnitude of this latter correction, while Figure ^ 
shows the corrected counts. 

Questions of sensitivity, completeness, spurious 
sources, and source confusion are important in a sur- 
vey exploited close to its limit, but are somewhat 
difficult to answer in view of the non-uniform noise 
properties of the mosaiccd images and the fact that 



some but not all of the sources are resolved. The ob- 
served field contains approximately 10^ independent 
beam areas and some 750 sources, resulting in over 
100 beam areas per source except in the few small re- 
gions of highest source density within the highest sen- 
sitivity area (only a few percent of the field). Hence 
we have not quite reached the sensitivity and source 
density at which confusion would become a serious 
problem. We have established that there are only a 
few deviations from a Gaussian distribution in the 
negative pixels, due to fluctuations similar to those 
edited out in SFIND. We thus expect approximately 
10 spurious sources due to noise statistics with our 
4ct peak flux cut-off. The catalogue completeness can 
be described in terms of the factors required to cor- 
rect for the primary beam weighting and the presence 
of resolved sources ( Hopkins ct al. 199^ ) . The survey 
has 50% completeness at that peak flux density which 
requires a factor of 2 correction for these effects. The 
peak flux density of the 50% completeness limit varies 
with radial position, r. It is 54 /LtJy for r < 12'; 69 /xJy 
for r < 21'; and 83 /iJy for r < 26'. The faintest cata- 
logued source has an integrated flux density of 45 /i Jy. 

3.2. Optical counterparts 

The optical catalogue is derived from Anglo- Australian 
Telescope (AAT) prime-focus CCD observations made 



in the Johnson-Kron-Cousins i?-band (Georgakakis 
|et al. 1999 ). The most probable optical candidate 
(where one exists) is chosen by searching a radius of 
5", and selecting the source with the least probabil- 
ity of being an accidental alignment (given the known 
surface density of sources as bright or brighter than 
the candidate) provided that that probability < 5%. 
While some true associations may be missed by this 
method it minimises spurious associations. Of the 
773 detected radio sources, 52% have been optically 
identified. As explained by Georgakakis et al. 1999 



the 50% completeness limit of the optical survey is 
R = 22.5 although sources fainter than this have been 
reliably detected and are included in the analysis. 

Figure ^ presents histograms of the distribution of 
apparent i?-band magnitude as a function of the radio 
source flux density. The median optical magnitude 
(^mcd) of the detected optical counterparts falls with 
decreasing radio flux density. The magnitude bound- 
ing the flrst decile (ii!io) also decreases to fainter opti- 
cal magnitudes as the radio brightness declines. The 
fraction of sources in each flux density bin with identi- 
fled optical counterparts, /, is smaller for the faintest 
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bin, to our fixed limiting optical magnitude. 
4. Discussion 

Over 50% of the 5*1.4 > 100 ^Jy radio sample is op- 
tically identified to i? = 22.5, with a lower rate (44%) 
for fainter radio sources. Visual inspection shows that 
the faint (< 100 /^Jy) sources are mostly identified 
with single optical galaxies, rather than interacting 
systems. This may reflect the fact that these identifi- 
cations are made quite close to the limit of the CCD 
images, or it may indicate that the very faint radio 
counterparts to a fixed optical magnitude are located 
in galaxies different from harbouring slightly brighter 
sub-mjy radio sources (c.f. Gruppioni et al. 1997). 

Figure |3| reveals considerable scatter in different 
determinations of the 1.4 GHz source counts below 
about 1 mjy. Counts from different surveys fluctuate 
by amounts large compared with the Poisson errors 
and with the respective image variances. The effect 
can be seen between counts in the 50' diameter survey 
area presented here and those in the entire 2° Phoenix 
field. It seems unlikely that limited dynamic range, 
sidelobe confusion, or the blending of faint sources to 
mimic single brighter sources could account for the 
effect within the Phoenix data set or within the other 
catalogues. 

It is known that large samples of radio sources 
brighter than w lOmJy at 1.4 GHz show no o r weak 
evidence of clustering (e.g. Baleisis ct al. 1998 ). How- 
ever, fainter sam ples may have reveal ed clustering 
( Cress et al. 1996 ), and it is predicted ( Benn fc Wall 
1995 ) that a 30' diameter survey with 5',„in < 100 ^Jy 



might well reveal fluctuations above Poisson statis- 
tics. Indeed, Benn & Wall caution that deep radio 
surveys with pencil beams may not be representative. 
Since the clustering scale responsible for the observed 
fluctuations is likely to be p ~ (100 Mpc)"^, ex- 
tending the Phoenix survey at its current sensitivity 
limit to a field of several degrees could provide con- 
straints to large-scale structure on scales intermediate 
between those probed by microwave background ex- 
periments and by optical redshift surveys. 

5. Conclusions 

A new 1.4 GHz radio survey with the ATCA has 
been completed as part of the Phoenix Deep Survey. 
A catalogue of 773 sources with 5*1.4 > 45 /iJy has 
been compiled. This sample is homogeneously se- 
lected and thus uniform across the entire flux den- 



sity range. The new sample provides an opportunity 
to investigate a large number of faint radio galaxies 
through multicolour photometry and spectroscopy. 

Statistical fluctuations in observed 1.4 GHz source 
counts below 500 /xJy provide a tantalising hint that 
very sensitive large area 1.4 GHz surveys may be 
used to provide constraints to cosmic structure on 
the scale p ^ (100 /i""'^ Mpc)^, intermediate between 
scales probed by microwave background experiments 
and by optical redshift surveys. 
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Fig. 1. — Grey scale image of the central 50' diameter 
region of the mosaic. The brightest source (in the NE 
of the image) is 23 mJy, the faintest 45 /iJy. The insert 
shows greater detail for a small portion of the image. 




Si4 (mJy) 

Fig. 3. — Normalised differential 1.4 GHz source 
counts. Crosses are a compilation of data from other 
surveys ( Windhorst et al. 1993| ). Open circles, open 
squares and filled triangles arc from radio observa- 
tions of the PDF. Filled s quares are from th e Hubble 
Deep Field South region ( Norris et al. 1999| ). 
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Fig. 2. — Points (right scale) show the raw source 
counts, with vertical error bars denoting Poisson vari- 
ation and horizontal bars denoting the range of source 
flux density averaged together. The fraction of the 
mosaic area covered by noise with 5a > Si a is shown 
as a dashed curve, and the weighting correction de- 
rived from this is shown a a solid curve. Approxi- 
mately 10% of the image has 5a < 50 ^Jy. 



Fig. 4. — Histograms of optical counterpart magni- 
tudes for selected ranges of 1.4 GHz flux density. The 
median i?-band magnitude of the optical counterparts 
in the given flux density range is shown for each panel. 
The fraction / of optical counterparts identified for 
the flux density range is also shown. The vertical 
dashed line indicates the completeness limit of the 
sample at i? = 22.5. 
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